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Determining the Dark Matter Relic Density in the Minimal Supergravity 
Stau-Neutralino Coannihilation Region at the Large Hadron Collider 
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We examine the stau-neutralino coannihilation (CA) mechanism of the early universe. We use the 
minimal supergravity (mSUGRA) model and show that from measurements at the Large Hadron 
Collider one can predict the dark matter relic density with an uncertainty of 6% with 30 fb _1 of data, 
which is comparable to the direct measurement by Wilkinson Microwave Anisotropy Probe. This 
is done by measuring four mSUGRA parameters mo, mi/2, Aq and tan/3 without requiring direct 
measurements of the top squark and bottom squark masses. We also provide precision measure- 
ments of the gaugino, squark, and lighter stau masses in this CA region without assuming gaugino 
universality. 

PACS numbers: 11.30.Pb, 12.60.Jv, 14.80.Ly 
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One of the important aspects of supcrsymmctry 
(SUSY), particularly when it is combined with super- 
gravity grand unification (SUGRA GUT) [2,0, is that it 
resolves a number of the problems inherent in the stan- 
dard model (SM). Aside from solving the gauge hierar- 
chy problem and predicting grand unification at the GUT 
scale M G ~ 10 16 GeV, subsequently verified at LEP 0, 
SUGRA GUT allows for the spontaneous breaking of 
SUGRA at the Mq scale in a hidden sector, leading to 
an array of soft breaking masses. The renormalization 
group equations then show that this breaking of SUGRA 
leads naturally to the breaking of SU(2) x U(l) of the 
SM at the electroweak scale, with SUSY breaking masses 
around a TeV for most of the SUSY parameter space. 

An additional feature of SUSY is that models with R- 
parity invariance give rise to a cold dark matter (CDM) 
candidate which is generally the lightest neutralino 
(X°). The CERN Large Hadron Collider (LHC) should 
be able to produce the and study its properties. Di- 
rect detection experiments for Milky Way DM would al- 
low for a determination of the DM mass and its nuclear 
cross section. If these are in agreement with the LHC de- 
termination of the Xi properties, it would help confirm 
the important point that the Milky Way DM was indeed 
the Xv However, this would not explicitly verify that 
the X? was the DM relic particle produced during the Big 
Bang. To do this, one would need to deduce the relic den- 
sity fl~oh 2 and compare with Ocdm'i 2 as measured as- 
Xi 

tronomically by Wilkinson Microwave Anisotropy Probe 
(WMAP) [|. 

In this Letter we describe a series of measurements in 
the stau-neutralino (ti-x?) coannihilation (CA) region 
where, in the early universe, the fi and the x? annihi- 
late together into SM particles, determining the relic DM 
abundance observed today. We show how to measure the 
sparticle masses, confirm we are in the CA region, mea- 
sure the SUSY parameters, and establish a prediction of 

fl~ h 2 . 
Xi 



To carry out this analysis it is necessary to assume a 
model that encompasses both LHC phenomena and early 
universe physics. Since the analysis is new and quite com- 
plicated we consider the simplest SUGRA model (mini- 
mal SUGRA or mSUGRA) [l[ with universal soft break- 
ing masses. However, we show below that it is possible 
to test experimentally gaugino universality; other non- 
universality models will be considered elsewhere. The 
mSUGRA model depends on one sign and four parame- 
ters: mo (universal sfermion mass) , m 1 / 2 (universal gaug- 
ino mass), Aq (universal soft breaking trilinear coupling 
constant) , tan (3 (the ratio of vacuum expectation values 
of two Higgs doublets), and the sign of /1 (the bilinear 
Higgs coupling constant). After we include all experi- 
mental constraints 0, the allowed mSUGRA parameter 
space with \i > (as preferred by b — > sj and the muon 
g — 2 [6]) has three distinct regions picked out by the 
CDM constraints 0: (i) the CA region where both m 
and TO1/2 can be small, (ii) the focus-point region where 
the x? bas a large Higgsino component and too is very 
large but m 1 / 2 is small, and (iii) the funnel region where 
both too and m 1 / 2 are large and the neutralinos can an- 
nihilate through heavy Higgs bosons (2M~o ~ Ma ,h q )- 

Al 

We consider here the CA region with fj, > 0. This re- 
gion is generic for a wide class of SUGRA GUT models 
(with or without gaugino universality). If the muon g — 2 
anomaly maintains, then the focus-point and funnel re- 
gions are essentially eliminated. 

The CA region has a striking characteristic of the fi 

AM : 



M fi - 



-M~o 
Xi 



-» TTXl 

c0 _ u 



; decays are 



and X? being nearly degenerate i.e 
(5-15) GeV. Thus, the x" -> rh 

dominant and the branching ratio for x 2 ~ * ~ y ^Xi is 
essentially zero (1= e or /1, and Ir is the lighter selectron 
or smuon). The existence of this near degeneracy would 
be a strong indication that we are in the CA region. 

In order to determine fl~oh 2 one must know all the 

mSUGRA parameters. In a previous study 0, TO1/2, too, 
Aq, and tan/3 were determined in the bulk region assum- 
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ing that it is possible to measure the gluino (g), squark 
(g), lighter bottom squark (sbottom or 61), Er, x®? an( i 
Xi masses in g — > 661 and — > (7X2 — * I^R ~^ Q^Xi 
decays, using "end-point" techniques Q. The determi- 
nation of Mg is very difficult if both g tti (here t\ is 
the lighter top squark or stop) and g — > 661 can occur [nj 
and the methodology of disentangling this background is 
not known yet. Also these techniques cannot be utilized 
for the CA case because the x° ~^ $R : - * ^Xi decay is 
essentially absent. 

While the CA region is particularly challenging, we 
show that it is indeed possible to determine all four pa- 
rameters accurately from measurements at the LHC. It 
has been recently shown 111, LL2J that the CA region can 



be established and that a measurement of AM can be 
made (provided the r identification can be done for visi- 
ble r p™ > 20 GeV) assuming A and tan/3 are known. 
The small AM value is experimentally characterized by 
a low energy r from a f\ — > TXi decay. With the ad- 
dition of some new datasets and variables, in particular 
with final state 6-quark jets, we show that we can (a) 
measure the g, q~L, X21 Xi; an d t\ masses in the case of 



M- q 



M 9 



> 



tion, 



without the mSUGRA assump- 
and (c) 



X2.X1 

(b) determine the mSUGRA parameters, 

predict f2~o/i 2 , which can be compared with the astro- 
Xi 

nomical determination of f^cDM^ 2 . Of particular note, 
our method effectively obviates the need to separate the 
final states arising from the third generation sparticles, 
such as stops (ii, £2), sbottoms (b\, 62), and staus (fi, 
T2). The procedure of extracting the model parameters 
is general and can be applied to other regions of the pa- 
rameter space or to more general SUGRA models. 

We select an mSUGRA reference point, shown in Ta- 
ble H where n~ h 2 = 0.10 and AM = 10.6 GeV. The 
Xi 

total production cross section at the LHC is 9.1 pb where 
the gq production has the largest contribution. Events 
are generated using ISA JET [13], followed by the PGS4 
detector simulation 



14] 



We analyze three samples with 
the final state of large transverse missing energy (fir) 
along with jets (j's), r's, and b's: (i) 2r + 2j + fir, (ii) 
ij + fir, and (iii) lb + 3j + fir- The kinematics in both 
2r + 2j + fir and the 16 + 3j + fir samples depend on 
all four mSUGRA parameters, while the 4j + fir sample 
is mostly sensitive to toq and m^. 



TABLE I: SUSY reference point (masses in GeV) 

350 GeV, m = 210 GeV, tan/3 = 40, A = 0, and fi > 0. 
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The primary SM backgrounds for the 2r + 2j + fir 
final state (and the other two samples) are from ti, 
W+jets, Z+jets and QCD production. The sample is se- 
lected using the following cuts [ll|: (a) N T > 2 (\r)\ < 2.5, 
p^ s > 20 GeV, but > 40 GeV for the leading r); (b) 
Nj > 2 (\ri\ < 2.5, E T > 100 GeV); (c) fir > 180 GeV 
and E T ]1 +E T j2 + fir > 600 GeV; and (d) veto the event 
if any of the two leading jets are identified as a b jet. In 
order to identify X2 ~ * — * TT Xi decays we categorize 
all pairs of r's into opposite-sign (OS) and like-sign (LS) 
combinations, and then use the OS minus LS (OS— LS) 
distributions to effectively reduce the SM events as well 
as the combinatoric SUSY backgrounds. We reconstruct 
the decay chains of qh — ► qx% — > qrfi — > qrrxi using the 
following five kinematic variables: (1) a, the slope of the 
distribution for the lower energy r in the OS— LS di-r 
pairs, (2) MP° ak , the peak position of the visible di-r in- 
variant mass distribution, (3) Mj^ 1 , the peak position of 

the invariant j-r-r mass distribution, and (4, 5) MP. eak , 
the peak position of the invariant j-r mass distribution 
where each r from the OS— LS di-r pair is examined sepa- 
rately. Note that we have used the peak positions instead 
of the end-points because of the r's in the final state. 

We follow the recommendation of Ref. for the 4j 
+ fir sample. The peak value, M^Z , of the variable 
M e ff = fir + J2 A j Cts Et j , which is a function of only 
the g and q masses, is reconstructed for each event that 
passes the following selection cuts: (a) Nj > 4 (|ry| < 2.5, 
E T > 50 GeV, but > 100 GeV for the leading jet); (b) 
fir > 100 GeV; (c) Transverse sphericity > 0.2; (d) Veto 
on all events containing an isolated electron or muon with 
p T > 15 GeV and \r]\ < 2.5; and (e) fir > 0.2M off . Again 
we require that none of these jets be identified as a 6 jet. 

Similar cuts are used to make the 16 + 3 j + fir sample. 
We introduce a new variable, M c ^ pcak , similar to M^ ak , 
but requiring that the leading jet be from a 6 quark. 

The measurement of a small value of a from the 2r + 
2j + fir sample indicates low energy r's in the final state 
(thus AM is small) and provides a smoking-gun signal for 
the CA region. In Fig. [TJ we show the distributions 
obtained by the OS— LS technique for various AM values. 
Note that a only depends on M^ and M^o (see Fig.[TJ. 

To get a set of measurements of the sparticle masses we 
use the remaining variables from the 2r + 2j + fir and 
4j + fir samples. The variables Mj^ and Mj 5 ^^ probe 
the q L -> qxl 



qr fi — > qrrxi decay chains. To help 
identify these chains we additionally require OS— LS di-r 

pairs with M TT < M™^ P om ^ and construct Mj TT for 
every jet with Et > 100 GeV in the event. With three 
jets, there are three masses: M^ T , M^ T , and M^ T , in 
decreasing order. We choose M^ T for this analysis 0]. 
Figure [2] shows the Mj TT distributions for two different 
qL masses, and Mj^ pcak as a function of Mg L and M^o, 
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FIG. 1: [top] The distribution of the lower-energy r's 

using the OS— LS technique in the three samples (arbitrary 

luminosity) of SUSY events with AM = 5.1, 10.6 and 16.9 

GeV, where only the fi mass is changed at our reference point. 

[bottom] The slope (as a in the text) as a function of 

the relative change of Mf 1 (therefore AM) or M-o from its 

Xi 

reference value with the other SUSY masses fixed. The bands 
correspond to statistical uncertainties with 10 fb _1 . 



keeping AM constant. Similarly, one can show that the 
^j2) P oak va j ue depends on the <jl, ^i an d Xi masses. 
The value of M p ^ ak , extracted from ij + frp sample, has 
been shown to be a function of only the ql and g masses. 

The determination of the sparticle masses is done by 
inverting the six functional relationships between the 
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FIG. 2: [top] The MfX distributions using the OS-LS tech- 
nique for SUSY events at our reference point, but with Mg L 
— 660 GeV (yellow or light gray histogram) and 840 GeV 
(green or dark gray histogram), where 748 GeV is our refer- 
ence point; [bottom] The peak position of the mass distribu- 
tion as a function of M^-o or M^ L . The bands correspond to 

statistical uncertainties with 10 fb" 1 . 



variables and the sparticle masses to simultaneously solve 
for the g, Xi ii Tl, and average qL masses and their un- 
certainties. The six parametrized functions are: AfP° ak 
= A (M-o, M-o, AM), a — / 2 (M- , AM), M^ peak 

A2 Al Al J 



^ M xV M # 

M-o, AM), and M e p ^ ak 



M (2)poak = U 6 (Mfj , M-o, 



Mg). With 10 



fb _1 of data, we obtain (in GeV) Mg = 831 ± 28, 

M~o = 260 ± 15, M~o = 141 ± 19, AM = 10.6 ± 2.0, 

X2 Xi—, 
and Mg L = 748 ± 25 [15J. The accurate determination 

of AM would also confirm that we are in the CA region. 

We also test the universality of the gaugino masses at 

the GUT scale. We measure M~JM~o = 5.9 ± 0.8 and 

A 1 

Mg/M ~o — 3.1±0.2, validating the universality relations 

to 14% and 6%, repectively. This non-trivial determi- 
nation of the additional gaugino masses along with the 
mSUGRA parameters require all six observables. The 
formalism developed here can work for other model with 
similar two-body decay processes. 

Since our primary goal is to determine 0,-oh 2 in the 

Xi 

mSUGRA model we next determine mo, m^, Aq and 

tan/3. M P g ak and Mj^ peak are insensitive to Ao and 
tan/3, and provide a direct handle on mo and mi/2 be- 
cause they depend only on the qi, (first two generations), 
g, x.2 aim Xi masses (see Fig. [3]). On the other hand, 
MP° ak and pcak provide a direct handle on Aq and 



tan/3. M r P r cak depends on M? t ; M o ( £ )pcak depends on 
Mjj and Afg , since both the t\ and b\ decays always 
produce at least one b jet in the final state. Figure [4] 
shows the values of MPf: ak and M 



(fr)peak 
off 



as functions of 



Ao and tan/3 since the off-diagonal elements of t\ and 
Si/fi mass matrices depend on M t (A t + jt/cot/3) and 
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FIG. 3: The dependence of Mj^ pcak (top) and M^ ak (bot- 
tom) as a function of mi/2 and mo. The bands correspond to 
statistical uncertainties with 10 fb _1 . 
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and without b jets. The accuracy of the SI -o/i 2 calcu- 

Xi 

lation at 30 fb _1 of data is expected to be comparable 
to that of Ocdm'* 2 by WMAP. This approach will allow 
us to determine the relic abundance at the LHC for any 
model where the CA is dominant in the early universe. 
Thus, it is possible to confirm that the DM we observe 
today were Xi's created in the early universe. 
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FIG. 4: The dependence of M ( £ )poak (top) and M T p T cak (bot- 
tom) as a function of tan/3 and Ao- The bands correspond to 
statistical uncertainties with 10 fb _1 . 



M b / T (A b / T + (j, tan/3), respectively. Combining these four 

measurements and inverting, we find mo = 210 ±4 GeV, 

mi /2 = 350±4 GeV, A Q = 0±16 GeV, and tan/3 = 40±1 

with 10 fb -1 of data pj|. 

After measuring the mSUGRA variables we calculate 

~oh 2 using DARKSUSY [16|. The calculation also involves 
Xi 

the Xi mixing matrix which we have determined in the 
mSUGRA case. In the CA region, Q,^ah 2 depends cru- 
cially on AM due to the Boltzmann suppression factor 
e -AM/k B T nl relic density formula 17]. Figure [5] 
shows contour plots of the la uncertainty in the fl^oh 2 - 

AM plane. The uncertainty on fl~oh 2 is 11 (4.8)% at 

Xi 

10 (50) fb _1 0. Note that it is 6.2% at 30 fb" 1 [Hj], 
comparable to that of the WMAP measurement [B| . 
In conclusion, we have established a technique for a 

precision measurement of tt-oh 2 at the LHC in the fi- 

Xi 

Xi CA region of the mSUGRA model. This is done using 
only the model parameters, determined by the kinemati- 
cal analyses of 3 samples of fk + j's (+ t's) events with 
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FIG. 5: Contour plot of the la uncertainty in the fL-o^i -AM 
plane with 10 and 50 fb _1 (outer and inner ellipses). 
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